Low-energy bound states at interfaces between superconducting and block 
antiferromagnet regions in K x Fe 2 _ySe 2 
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The high-T c alkali doped iron selenide superconductors YL x ¥e2- y Se2 have been recently shown to 
be intrinsically phase separated into Fe vacancy ordered block antiferromagnetic regions and super- 
conducting regions at low temperatures. In this work, we use a microscopic five orbital Hubbard 
model to obtain the electronic low-energy states near the interfaces between block antiferromagnets 
and superconductors. It is found that abundant low-energy in-gap bound states exist near such in- 
terfaces irrespective of whether the superconductor has d- or s-wave pairing symmetry. By contrast, 
it is shown how nonmagnetic scattering planes can provide a natural means to distinguish between 
these two leading pairing instabilities of the K a! Fe2- ! /Se2 materials. 
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I. INTRODUCTION 

Alkali doped iron selenide materials 1 -^ A x ¥e2- y ^>e2 
(A=alkali element) are high temperature superconduc- 
tors with T c ~ 32K— that also undergo a transition to a 
iron vacancy ordered structure at Ts ~ 578K—~— followed 
closely by a transition to a block anti-fcrromagnctic state 
(BAFM) with Neel temperatures T N ~ 559K£rl The 
BAFM develops within the -\/5 x y/5 arrangement of iron 
vacancies with a magnetic moment of m ~ 3.31/^s/ Fe. 
This is the highest reported magnetic moment among 
all iron pnictide and iron chalcogenide superconductors. 
The existence of the above phases within the super- 
conducting state has generated intense debate about 
the interplay between these seemingly competing phases, 
and about the origin of superconducting pairing leading 
to high temperature superconductivity in this unusual 
environment £ 

The presence of phase separation between supercon- 
ducting regions and the BAFM has been confirmed us- 
ing a variety of experiments including spectroscopy* 9 - - — 
microscopy* 9 *!^— and x-ray diffractio n 16 ' 18 ' 19 tech- 
niques. The phase separation has been observed along 
the c-axis and also in the FeSe plane as nanoscale fil- 
amentary superconducting domains existing within the 
vacancy ordered BAFM regions! 11 ' 12 ' 18 Though the mor- 
phology and physical properties of the metallic nano 
domains can depend on the thermal histor y 16 ' 19 and 
stoichiometryfi^ a number of alkali iron selenide samples 
have been observed to be free from iron vacancies in the 
superconducting regions! 10 ' 15 ' 20 

Angular resolved photocmission (ARPES) experiments 
investigating the metallic regions observe a Fermi surface 
with two electron pockets around the M points and a 
small electron pocket developing around the Z pointi 2 ^- 2 ^ 
The absence of a hole pocket around the T point has rein- 
vigorated the discussion of the origin of superconductiv- 
ity in the Fc-based materials. Specifically, the absence 



of (tt, 0) nesting between hole and electron Fermi surface 
sheets prevents the leading s± pairing instability found 
e.g. within spin fluctuation mediated pairing^ Theoreti- 
cal calculations for the superconducting ground state are 
currently inconclusive and predictions have been made 
for both rf-wav o 25 ' 26 and s-wav o 27 ' 28 pairing symmetry. 

Measurements within the superconducting state find 
a nodclcss superconducting gap i 14 ' 21 ' 22 ' 29 Although 
ARPES detects a nearly isotropic gap structure* 21 - STM 
measurements have observed a double gap feature in the 
local density of states (LDOS)i 13 ' 14 This does not iden- 
tify the gap symmetry in alkali iron selenides since nei- 
ther d-wave nor s-wave symmetry possess any symmetry 
enforced nodes at k z = for this Fermi surface. How- 
ever, ARPES measurements in Ref. |30| find an isotropic 
superconducting gap structure at the Z point where a 
electron pocket exists around the T point indicating an 
s-wave gap symmetry. 

From the above discussing, it is evident that new the- 
oretical methods to determine the pairing symmetry of 
these materials are desirable. In this paper we study 
the in-plane interface states near boundary regions of a 
BAFM with ordered arrangement of iron vacancies and 
a next nearest neighbor (nnn) s-wave or nearest neigh- 
bor (nn) paired d-wave superconductor. Note that in 
the absence of a Fermi pocket around the T point the 
nnn s± state does not have a sign change between Fermi 
pockets. We denote this superconducting state by s-wave 
with the understanding that it represents nnn supercon- 
ducting pairing interaction in real space. We use a five 
orbital tight binding Hamiltonian and include supercon- 
ductivity by generating pairing from a spin fluctuation 
exchange mechanism within an RPA weak coupling the- 
ory. The model is able to reproduce the electronic struc- 
ture of K a; Fe2-j / Se2 and stabilise a BAFM structure as 
well as a superconducting state with s-wave or d-wave 
symmetry. It is found that the presence of an interface 
between these two ordered phases leads to the forma- 
tion of in-gap bound states which should be detectable 



by future STM measurements. Bound interface states 
are generated for both s-wave and <i-wave symmetry due 
to the magnetic nature of the BAFM. A qualitative dif- 
ference between s- and d-wave order is shown to arise 
near e.g. (110) non-magnetic scattering boundaries such 
as cracks, grain boundaries or free surfaces within the 
superconducting phase. 

The paper is organized as follows, in the next section 
we provide details of the model used for our calculations. 
This is followed by an analysis of the electronic structure, 
magnetic order, bulk superconducting order and finally 
we discuss the results for the LDOS in the simplified in- 
homogeneous situation of an interface between a BAFM 
and a superconductor. 
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II. MODEL 

The five-orbital model Hamiltonian is given by 

H = Ho + Hint + HsC + H vacancy- (1) 

The first term is determined by the tight-binding band 
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Here the operators c\ {c\^ a ) create (annihilate) an elec- 
tron at the i-th site in the orbital /i and with spin pro- 
jection er, and /zo is the chemical potential. The indices 
/i and v correspond to the d xz , d yz , d x 2_ y 2, d xy and 
d 3z 2 iron orbitals. The hopping parameters t^ are ob- 
tained by slightly modifying the tight-binding parame- 
ters of KFe2Se2 in Ref. [3l|. This is similar to the pre- 
scription used by Maier et al£^ and consists of shifting 
the d xz /dy Z orbital nearest neighbor hopping by -0.08 
cV and d xy orbital hoppings by 0.08 cV. The resulting 
band structure and Fermi surface are shown in Fig. [T] 
It can be seen by comparing with the KFe2Se2 bands in 
Ref. [l9| that the modification of the hopping parameters 
shifts down the hole band near the T point and leads 
to a Fermi surface with only electron sheets. The Fermi 
surface shown in Fig. [TJb,c) corresponds to an electron 
doping of K x Fc2-j,Se2 with x = 0.3, and y — (0.15 elec- 
trons per Fe) which may be contrasted to the undoped 
compound with x = 1.0, and y = (0.5 electrons per 
Fe). The Fermi surface calculated in Fig. [TJc) contains 
a small electron pocket around the Z point in addition 
to the Fermi pockets around the M point. The LDOS 
calculated at k z = in Figfljd) contains mainly d xz , d yz 
and d xy orbitals near the fermi energy whereas the con- 
tribution from the d x 2_ y 2 and d 2 z orbitals near the Fermi 
energy arc negligible. 

The second term of Eq.fl]) describes the intra-site 
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FIG. 1. (Color online) (a) Band structure of K;cFe2-;;Se2. 
Here the color codes represent green: d xz , red: d yz , cyan: 
d x 2_ y 2, blue: d xy , yellow: d z 2. (b) Fermi surface calculated 
at k z = 0. (c) Fermi surface at k z — it. (d) Orbital resolved 
LDOS at k z — in the normal state at an electron doping of 
n = 0.15 e~/Fe. 



Coulomb interaction 



Hint = U 2J ni^riiui 



{u 1 - — ) 2_j n ^ n ™<?' 



i,/iO,erer 



(3) 



2J£ S^-S^ + J* J2 i^Ci 



I/cT L -ll/cTi 



which includes the intra-orbital (inter-orbital) interaction 



U ill'), the Hund's rule coupling J and the pair hopping 
energy J'. We will assume orbital and spin rotational 
invariance where the relations U' = U — 2 J and J' = J 
hold. 

The third term in Eq. ([T]) is the superconducting pairing 
term 
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ing phase could be stabilised by additional pair scattering 
processes than the ones allowed by Fig. QJb) caused by 
e.g. grazing hole pockets near the T point, or important 
inter-band processes between hole and electron processes 
at k z = 7Tii 

The last term in the Hamiltonian ([1} incorporates the 
iron vacancies. We model the vacancies as impurities 
with sufficiently large onsite potentials that remove any 
electronic states from their location 



with the order parameter Ajt" = r^'^(i, j) (cj^Ci^-j-). 
Here, r^(i,j) is the strength of an effective attraction 
that is generated from a spin fluctuation exchange mech- 
anism. To this end, we follow the procedure given e.g. 
by Graser et aZ.— and calculate the singlet pairing vertex 

i*Kk,*>, u ) = [\u s x f PA {k - fc', w )C/ s + \w 

-\u c X$ PA {k-k>,w)U°+\u%l. (5) 

Note that we have retained only those pairing vertices 
of the form r££(i, j). This is done only for computational 
speed since we do not find any qualitative changes in our 
results by including additional inter-orbital interactions 
as they are quite small in magnitude. The pairing inter- 
action strengths are fixed by requiring that U and J gen- 
erate a superconducting gap magnitude of A = 12mcV 
for both the (i-wave and s-wave symmetry which agrees 
reasonably well with gap magnitude observed e.g. in 
ARPES experiments^ The obtained pairing symmetry 
is determined simply by the Fermi surface used as input 
to Eq.©, i.e. the Fermi surface in FigQJb) [FiglHc)] 
generates d [s] pairing, respectively. In Table U we show 
the self-consistent gap magnitudes A^ on each orbital. 
It can be seen that the gap primarily resides on the d xy 
orbital which we attribute to the intra orbital pairing in- 
teraction of the d xy orbital (r^|^|(r, r')) being at least 
twice that of the d xz /d yz orbital (T% z z >% z z {y, r')). 

We stress that the present calculation is 2D, and for 
the results discussed below we use a tight-binding band 
with the Fermi surface shown in Fig. [TJb). It is only 
for generation of the pairing couplings T that we utilise 
the fact that d- or s-wave pairing is preferred within the 
RPA spin-fluctuation approach by Fermi surfaces corre- 
sponding to Fig.QJb) and Fig. HJc), respectively. This is 
simply a way to simulate an s-wave state in the presence 
of the Fermi surface of Fig. [TJb). Such a superconduct- 
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which adds the potential V vacancv at the vacancy sites i* . 
After a standard mean-field decoupling of the onsite 
interaction term (|3]) we arrive at the following multiband 
Bogoliubov de-Gennes equations^! 
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The local densities and the SC order parameters are ob- 
tained sclf-consistcntly through iteration of 
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In the following section, we discuss the results of the 
above procedure for the magnetic, and superconducting 
properties applicable to iron chalcogenide superconduc- 
tor K r Fe2-ySe 2 . 



III. RESULTS 
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(1,0), (-1,0) 


0.62 


0.69 


0.03 


7.19 





d-wave 


(0,1), (0,-1) 


-0.69 


-0.62 


-0.03 


-7.19 





s-wave 


(1,1), (-1,-1) 


0.68 


0.68 


0.02 


8.68 





s-wave 


(-1,1), (1,-1) 


0.68 


0.68 


0.02 


8.68 






TABLE I. Real space self-consistent gap values AJ 1 ." on each 
orbital for nearest neighbor and next nearest neighbor pair- 
ings. Units are in meV. 



The band structure properties of the iron based super- 
conductors generally include contributions from all five d- 
orbitals near the Fermi surface. A simple electron count 
and ARPES results reveals that the 122 iron selenides 
are strongly electron doped. We impose a 15% electron 
doping on the metallic region by adjusting the chemical 
potential to fi = — 0.23eV. The evaluated Fermi surface 
shape for k z — and k z = it are shown in Fig. [TJb) and 
0Jc), respectively. They agree well with ARPES observa- 
tions which find only electron pockets at the M points for 




:■■■■ 



FIG. 2. (Color online) (a) Orbital resolved LDOS of the 
BAFM region with a normal state band corresponding to 
Fig- E (b) Real-space map of the magnetic moment in the 
BAFM state. White: vacancy sites, green: 3.63/is/Fe, blue: 
-3.63/^s/Fe. (c) Same as (b) but for orbital order n xz -n yz . 
Blue: -0.011, green: 0.011. 



k. = and a small additional electron pocket around the 
Z point 3ir— In our calculations for the normal metal- 
lic state, the onsite Coulomb term has been fixed by 
U = 0.6eV and J = 0.25*7. 

In Fig. []Jd) we snow the orbitally resolved LDOS. It 
can be seen that the primary contributions to the den- 
sity of states near the Fermi level are from the d xy and 
dxz/dy Z orbitals. Though the orbital content can depend 
on the doping level as well as the presence of iron vacan- 
cies, observations on similar alkali doped iron selenide 
systems identify the above \>i g orbitals to be the domi- 
nant contributors of the Fermi surface at k z = 0.— 
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FIG. 3. (Color online) (a) LDOS in d-wave superconducting 
state. Black: total LDOS, blue: dxy, red: d xz /d yz . The inset 
shows the gap anisotropy around the Fermi surface centered 
at (0, 7r). The ir-axis is angle measured with respect to the k x 
direction, (b) Same as (a) but for the s- wave symmetry. 



discussed above. An electron density of 6 e~/Fe is ob- 
tained by a Coulomb interaction strength of U = 1.35eV 
and J = 0.25*7 in the BAFM state. The stoichiometry 
for the BAFM with V5 x y/E can be considered to be 
K0.8Fe1.6Se2 from an electron and vacancy concentration 
count. It has been previously argued by Luo el al£^ 
that a correct description of the magnetic state requires 
a larger Coulomb interaction term compared to what is 
typically used for pnictides. In addition, we note that 
in the phase separated region the normal metal is likely 
to have a smaller effective on-site Coulomb repulsion as 
otherwise the normal region would undergo a SDW tran- 
sition to a (jr, 7r) magnetic state and no such magnetism 
has been observed. 

Using the self-consistency approach discussed in sec- 
tion 2 we find a stable BAFM state with a magnetic mo- 
ment of 3.63 /is/Fe. The magnetic moment obtained is 
shown in a real-space map in Fig. I^b) . Calculation of the 
associated LDOS yields an insulating gap of around 0.35 
eV as shown in Fig.(2ja). This gap agrees reasonably well 
with the gap magnitude obtained in STM measurements 
of the BAFM regions^ Another result that comes out 
naturally from this calculation is the presence of orbital 
order defined as n xz — n yz as shown in Fig. [2jc). The 
origin of the orbital order is the orbitally selective sup- 
pression of hopping of an electron in a d xz (d yz ) orbital 
in the horizontal (vertical) direction next to a vacancy 
location.— 



B. Block antiferromagnetic state 



C. Superconducting state 



The presence of an ordered arrangement of iron vacan- 
cies with a v5 x y5 arrangement in the lattice struc- 
ture of Ka;Fe2-ySe2 stabilises a block antiferromagnetic 
state (BAFM). STM measurements, for example, find 
that these regions are insulating^ Additionally, experi- 
ments find an electron occupation of 6e~/Fei 15 ' 20 

In our model calculations, the chemical potential for 
the BAFM state is the same as in the normal metal 



As described in the model section, superconducting 
order is stabilised by including the pairing interactions 
generated within a standard weak coupling RPA method. 
Before turning to the discussion of the LDOS near bound- 
ary regions between superconducting and BAFM regions, 
it is instructive to briefly display the properties of the 
homogeneous superconductor. In Fig. [3] we show the 
orbitally resolved LDOS for the superconducting state 
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(b) Orbital Order 
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(c) Superconducting Gap 



|A| : 





35 40 45 50 55 60 



FIG. 4. (Color online) Spatial dependence of the order pa- 
rameters near an interface between a BAFM and d-wave su- 
perconductor. The system size is 80 x 10 with a BAFM of 
10 x 10 sites placed between site x — 40 and x = 50. (a) Map 
of the magnetic moment (/is/Fe) shown between sites x = 30 
to x = 60 and y — 1 to y — 5. (b) Same as (a) but for the 
orbital order (n xz — n yz ). (c) Orbital resolved averaged su- 
perconducting gap magnitude shown between x — 30, y — 10 
toi = 60, y — 10. The color scheme is the same as in Fig. [T] 



with d-wave and s-wave pairing symmetry, respectively. 
Due to the same gap amplitude in the two cases, and 
the fact that the LDOS is not phase sensitive, the final 
LDOS of the homogeneous phase is very similar in the 
two states as expected. The calculated LDOS are typi- 
cal for a nodeless superconducting state and display an 
apparent "two-gap" behavior similar to what has been 
observed in recent experiments ! 13 ' 14 The two sets of co- 
herence peaks close to ~ 12meV and ~ 25meV in Fig. [3] 
results from a significant gap anisotropy on the electron 
pockets as seen from the insets in Figs. |3] The super- 
conducting pairing interaction has been chosen to result 
in minimum gap amplitudes of about 12mcV which is 
of the order of experimentally observed amplitudes that 
have been observed to be around lOmcV. 



D. Interface between BAFM and superconductor 

1. BAFM/SC interface 

As mentioned in the introduction, the alkali doped iron 
selenide materials A 2: Fe2-, / Se2 are phase separated, and 
interfaces between BAFM and superconducting regions 
are abundant and an important situation to understand 
in detail. For such an interface, it can be seen from 
Fig. B|a) that the magnetism penetrates only a few lat- 
tice spacings into the superconducting region due to the 
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FIG. 5. (Color online) (a) Evolution of the LDOS upon ap- 
proaching a BAFM interface from x = 1 to x = 42 for s-wave 
gap symmetry, (b) Same as (a) but for d-wave symmetry of 
the superconducting order parameter. 



significant amplitude of the magnetisation in the BAFM 
regions. Since the BAFM is insulating, the superconduct- 
ing order parameter also does not exhibit any significant 
leaking into the BAFM regions as seen from Fig. B|c). 
This short-range proximity effect is qualitatively consis- 
tent with the existence of filamentary superconducting 
channels near these magnetic regionsj 11 ' 12 ' 18 The orbital 
order near the interface is shown in Fig. HJb). As seen, 
at particular sites near the interface, the orbital order 
undergoes significant enhancement by a factor of about 
ten compared to the size of the orbital order in the bulk 
BAFM. This takes place adjacent to the vacancy sites 
bordering the superconductor causing significant local 
LDOS modulations at these sites. Note that the abso- 
lute change in the orbital order near the interface is still 
weak with n xz — n yz ~ 0.11. 

What are the low-energy states near the BAFM and 
superconductor interfaces? In Fig. [5] we show the evolu- 
tion of the LDOS upon approaching the interface for both 
the s- and d-wave cases. As seen in both cases, prevalent 
low-energy in-gap bound states clearly exist, and there 
is no qualitative difference in the LDOS caused by the 
pairing symmetry of the bulk superconductor. This is 
contrary to the case of a single nonmagnetic impurity 
where Zhu et alM- showed that the d-wave state leads to 
in-gap bound states as opp osed to the nnn s-wave state 
(see however Ref. U\ and [43 for a more general discussion). 
The magnetisation of the BAFM is crucial for the genera- 
tion of in-gap states in both Fig. (5{a) and Fig. [SJb) since 
it breaks time-reversal symmetry and generates bound 
states for both s- and d-wave superconductors. Certainly, 
in a picture where one thinks of the interface as a row of 
magnetic and non-magnetic impurities the resulting in- 



(a) s-wave 



(b) d-wave 




grain boundary leads to the generation of a zero-energy 
resonant state due to a sign change of the gap function 
between the incoming and scattered momentum vector 
within a semiclassical picture4£~— In the present case, 
a similar scattering geometry which probes the different 
gap signs of the electron pockets should provide a qual- 
itatively different low-energy LDOS between the s- and 
d-wave superconductors. As shown in Fig. [51 a (110) 
wall placed at 45 degrees to the x-axis indeed leads to 
a bound state solutions for the d-wave state unlike the 
case for the s-wave gap symmetry. The origin of this 
difference is exactly the same as for the standard one- 
band d-wave materials 4&~— Thus, using scanning probes 
to locate such dislocations in the FeSe planes, local STM 
measurements near such regions could provide important 
clues to the pairing symmetry of these materials. 
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FIG. 6. (Color online) An 80 by 10 system with an interface 
between an s-wave nnn superconductor and a BAFM with 
vacancies forming a (110) wall. Here, U = 0.6, and J = 0.25U 
which leads to a nonmagnetic interface, (a) LDOS plots for 
x — 1 to x — 40 for s-wave gap symmetry, (b) Same as (a) 
but for d-wave gap symmetry, (c) Incident wave vector k 
at the (110) interface and reflected wave vector k' for s-wave 
symmetry, (d) Same as in (c) but for d-wave symmetry. Blue: 
IAI, Red: -IAI. 



IV. CONCLUSIONS 

We have provided a microscopic real-space study of 
the interplay between BAFM and superconducting do- 
mains relevant to alkali doped iron selenide materials 
Aa;Fe2- y Se2 which have been shown to exhibit signifi- 
cant in-plane and out-of-plane phase separation. Our 
modelling includes all five iron d-orbitals and the on- 
site Coulomb interaction treated at the mean field level. 
Superconductivity is incorporated through the spin- 
fluctuation RPA approach. A main result is that low- 
energy in-gap bound states should exist near interfaces 
between BAFM and superconducting regions irrespective 
of whether the pairing symmetry of the superconductor 
is d- or s-wave. By contrast, a non-magnetic (110) scat- 
tering plane only causes low-energy bound states for the 
d-wave pairing and should provide a means to distinguish 
between d and s-wave pairing symmetry of these materi- 
als. 



tcrference of the multiple generated bound states should 
generate an in-gap band of states leading to LDOS results 
similar to those shown in Fig. [5] 44 ' 45 



2. Scattering off a (110) wall 

In the case of cuprate d-wave superconductors it is 
well-known that a (110) diagonal surface/interface or 
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